Ferromagnetic spin-polaron on complex lattices 
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We present a simpler derivation of the exact solution of a spin-polaron in a ferromagnet and 
generalize it to complex lattices and/or longer range exchange interactions. As a specific example, 
we analyze a two-dimensional Mn02-like lattice (as in the ferromagnetic layers in LaMnOa) and 
discuss the properties of the resulting spin-polaron in various regimes. At strong couplings the 
solution is reminiscent of the Zhang-Rice singlet, however the electronic wavefunction involved in 
the singlet is dependent on the momentum of the singlet, and multiple bands may appear. 

PACS numbers: 71.10.Fd, 71.27.+a, 75.50.Dd 



I. INTRODUCTION 

The motion of a charged particle in a magnetically or- 
dered background is one of the central basic problems 
encountered in doped, magnetically ordered insulators. 
Well known examples are hole or electron doped ferro- 
magnetic insulators such as EuO,^ hole doped parent 
compounds of the colossal magneto-resistance materials,^ 
hole and electron doped parent compounds of high- 
temperature superconductors r'i^ etc. 

The exact solution of the problem in a two-dimensional 
(2D) antiferromagnetic lattice as represented by the 
cuprates still eludes us. Diagrammatic Quantum Monte 
Carlo calculations on an assumed Neel-ordered lattice 
and in a single band i J-like model provide an exact nu- 
merical solution to this approximated systemi^ In the 
real system, however, the hole propagates in an O 2p 
band while the spins are a result of a half-filled Cu (1^2 _y2 
band with a large Hubbard U. 

The other examples like EuO or LaMnOa are either 
ferromagnetic or, as in the case of LaMnOa, have ferro- 
magnetic 2D layers. As we show here, for these systems 
there is an exact solution available for a hole or electron 
propagating in either the same, or a different band from 
that of the spin background, and for any sign and magni- 
tude of the coupling. Exact solutions of this kind provide 
important information on the existence range of bound 
spin-polaron states as envisioned for example for Zhang- 
Rice (ZR) singletS)^ in which the hole propagates in an 
O sublattice and the local spins are on a transition metal 
(TM) sublattice. It is important to note that aside from 
cuprates, there is evidence that doped holes in mangan- 
ites LaMnOs and cobaltates Naa;Co02 also propagate on 
the O sublattice but are strongly coupled to the TM. 

Indeed, it has been known for a long time^i^i^ i^^i^^i^^ 
that an exact solution can be found for the Green's func- 
tion of a particle (electron or hole) moving in a lattice of 
ferromagnetically (FM) ordered spins, at zero tempera- 
ture. If the spin of the particle is parallel to the FM order, 
the solution is trivial: its energy is simply shifted by its 
exchange coupling to the FM spins. However, if the parti- 
cle spin is antiparallcl to the FM order, it can scatter and 
spin-flip by creating a magnon. Depending on the values 
of the various parameters and the total energy, this can 



result in a finite lifetime (incoherent scattering leading to 
a broad spectral weight) or in an infinitely-lived quasi- 
particle - the spin-polaron - comprising the bound-state 
of the particle and the magnon. In the later case, the 
dispersion of the spin-polaron can be significantly renor- 
malized compared to that of the free particle. 

Similar physics is found in a seemingly more com- 
plex Kondo-Anderson model that also admits an exact 
solutionJ^ This is not surprising, because the restriction 
placed there on the allowed occupation numbers of the 
d orbitals essentially maps it back to an electron inter- 
acting with FM-ordered spins. Effort has also been fo- 
cused on trying to extend this type of solution to finite- 
temperature (i.e. presence of multiple magnons), to finite 
particle concentrations, etc. While, to our knowledge, no 
other exact solutions have been found, such work has re- 
sulted in various approximations for the self-energyi^^ 

One common aspect of all these exact solutions are 
the assumptions (1) that the particle moves on the same 
lattice that hosts the FM spins (usually a simple cubic- 
like lattice in d dimensions, although generalizations to 
other cases are trivial), and (2) that the particle-spin 
exchange is local, i.e. purely on-site. 

In this article we show that the T = exact solution 
can be generalized to systems where these restrictions 
are lifted. In other words, to cases where the particle 
moves on the same lattice that hosts the FM spins but 
the exchange is longer range, as well as to cases where 
the particle moves on a different sublattice than the one 
hosting the FM spins so that the exchange is necessarily 
not on-site. We focus on a specific problem of the later 
type, and briefly comment on other possible generaliza- 
tions later on. While our formalism is similar in spirit to 
that used in Refs. [tUsHqI IToIITiIIT^ . it is in practice much 
simpler to use and more transparent. This is essential to 
allow us to find these generalized exact solutions. 

Such calculations are necessary in order to under- 
stand quantitatively direct /inverse angle-resolved photo- 
emission spectroscopy (ARPES) (when the particle is a 
hole/electron), in insulators which order ferromagneti- 
cally, for instance oxides like EuO, CuCr2S4, CuCr2Se4, 
MgFe204, etc. They may also be relevant to some extent 
for itinerant ferromagnets, given that spin-polarized elec- 
tron loss spectroscopy on thin Co films^^ reveals good fits 
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FIG. 1: Sketch of the lattice and the indexes of sites in the 
unit cell. The bottom-right cell shows the signs of the p orbital 
lobes, which define the signs of the hopping integrals ii,t2. 

of the measured spin-wave spectrum to effective Heisen- 
berg models, however more work needs to be done to 
understand how to properly extend it in this direction. 

Our results suggest possible spintronic uses for these 
materials, in terms of transport of spin-polarized cur- 
rents, since we prove that charge carriers with spin- 
polarization antiparallel to that of the FM background 
can propagate coherently, in other words scattering on 
magnons is not necessarily leading to a finite lifetime. 
Equally importantly, the results also give some useful in- 
sights on ways to improve our understanding of the prop- 
agation of electrons or holes in antiferromagnetic back- 
grounds. These issues are discussed in more detail below. 

The paper is organized as follows: in section II we in- 
troduce the specific model to be solved. In section III we 
give its exact solution, and comment on various possible 
generalizations. In section IV we present and analyze a 
selection of interesting results. Section V contains our 
summary and conclusions. 



II. MODEL 

We study the 2D model of Fig. [1] It consists of two 
sublattices: the one hosting the FM-ordered spins is a 
simple square lattice of lattice constant a, whose sites 
are marked by black squares (the analogs of the Mn sites 
in a Mn02 lattice). The charge carrier is moving on the 
other sublattice, which includes all the sites marked by 
circles (the analogs of O sites in a Mn02 lattice). The 
spin of the charge carrier is coupled through exchange to 
the spins on its two nearest neighbor sites. 

Such a model would describe for example a Mn02 layer 
of the parent compound LaMnOa if it was modeled as a 
charge-transfer insulator^iS, The Mn ions are in the id^ 
{S — 2) configuration, and have ferromagnetic order in 
each layer. In a charge-transfer insulator, doping would 
introduce a hole in the O 2p orbitals instead of emptying 
an Cg Mn 3d orbital. The spin of the hole would be 
AFM coupled to the neighboring Mn spins, with a Jq ~ 



|tp_dP/Ai5, where tp-d would measure the hybridization 
of the orbitals, and AE would be the overall energetic 
cost to move the hole to a Mn site. 

This model is also reminiscent of the Cu02 layers of 
cuprate parent compounds, which are charge transfer ma- 
terials with holes going into O 2p orbitals. In this case, 
S = 1/2 for the configuration of the Cu in the insu- 
lator, however of course these spins order antiferromag- 
netically, not FM as assumed to be the case here. 

In any event, our primary motivation here is to exem- 
plify the general solution for these types of problems, and 
we chose this model Hamiltonian because it is complex 
enough to demonstrate the full power of our solution, 
yet simple enough not to make the notations too cum- 
bersome. It also unveils some very interesting physics. 
The solution can be directly generalized to a much wider 
class of similar problems, as discussed below, and can 
therefore be used to describe realistic systems. 

We assume a total of N unit cells, with periodic bound- 
ary conditions. In all our results we let N — > oo. Each 
unit cell is indexed by a pair of integers and con- 

tains three sites: a spin site located at Rij = iax + jay 
and 2 incquivalent electron sites, denoted 1 (for the site 
on the x-rung, location Rij.i — Rij + §i) and 2 (site on 

the y-rung, location Rij,2 — Rij + f y). 

The hopping between the various electron sites is also 
indicated in Fig. [T] For simplicity we limit ourselves to 
nearest-neighbor hopping - generalizations are straight- 
forward. If the orbitals occupied by the charge carrier are 
p-orbitals (as drawn), then one must take ti = —t2 = t, 
whereas if they are s-type orbitals (not shown), one 
should choose ti — t2 — t. In any event, we will use 
general ti and <2 values for the derivation. The differ- 
ence between a hole vs. an electron charge carrier is the 
sign of t: t > for electrons, and t < for holes. At first 
sight one would think that the sign of t is irrelevant, since 
the model is particle-hole symmetric. In fact, while the 
energetics is the same, the sign is important for the wave- 
functions and therefore has interesting consequences. 

In the following, we assume that the charge carrier is 
an electron, and we will comment on the differences for 
holes where appropriate. 

We introduce the Fourier-transformed operators 

where c]^- ^ ^ creates an electron with spin a at site A = 
1,2 of the unit cell. In terms of these, the nearest- 
neighbor hopping of the electron is written simply as: 

^-E^.-(ii,.-M,. + H (2) 

where the sum is over the allowed values of k in the first 
BriUouin zone (BZ) (-f , f ] x (-f , ^], and 

e,^ = -2h cos _ 2t2 cos fctM^ . (3) 
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The hopping Hamiltonian is diagonalized trivially by us- 
ing the new operators 



1 / 



in terms of which the hopping Hamiltonian is: 



T : 



(4) 



(5) 



k,{T 



As expected for this two-site unit cell, there are two 
bands for the free electron (they happen to touch, there 
is no gap between them). If ti ~ t2 (s-orbitals), the 
ground-state is at the F-point = ky = 0. If ti — ~t2 
(p-orbitals) , the ground-state is moved to k^ = ky — n / a 
or equivalent points at the corners of the Brillouin zone. 
In fact, the whole dispersion is just translated by ^(1, 1) 
inside the Brillouin zone, hence the resulting physics is 
exactly the same in both cases. Differences between s- 
and p-orbitals (besides this overall shift) are only appar- 
ent if there is longer range hopping. Since we do not 
consider longer-range hopping here, in the following we 
will report results for the s-orbital case ti = t2- 

To the hopping term, we add the FM Heisenberg ex- 
change between the spins of magnitude S : 



n.. 



spins 



(6) 



where the sum runs over all units cells (we limit our- 
selves to nearest-neighbor exchange since generalizations 
to longer-range exchange are trivial, so long as there is 
no frustration) and the electron-spin exchange: 



+ Sii,2 • ( Sij + Sij + ] 



(7) 



where Sy,A = J2aj 4j,\,a^'4j,\,f3 ^^e spin of the elec- 
tron at site A = 1, 2 of the {ij) unit cell. 

The total Hamiltonian is the sum of the three terms of 
Eqs. (O, ([6]) and ([7]). For later convenience, we divide it 
in two parts: 



Ti = Tin 



where 



TYo — T - J ^ [Slj Sl^ 



V 



QZ 



2S' 



(8) 



(9) 



h3 



includes the hopping and the diagonal {zz) part of the 
FM Heisenberg interaction between the spins. V includes 
the remaining terms, i.e. the xy part of the FM exchange 
between spins, and the coupling of the electron spin to 
the spins located on its two neighboring sites. 



III. THE GREEN'S FUNCTION OF THE 
SPIN-POLARON 

We consider the case where a single electron is in the 
system and T = 0. Let |FM) = \ + S,+S,..., +S) be the 
ground-state of the spins, in the absence of the electron. 
If the electron has spin-up the problem is trivially solved, 
since the xy parts of all exchanges have zero action in 
this subspace (no spins can be flipped). We therefore 
only treat explicitly the case of a spin-down electron. 

We define 



A,',( 



(10) 



and introduce the 2x2 Green's function matrix: 

G{k,Lo) = (FM|vI/^,^G(^)vI't JFM), (11) 

where the resolvent is G{uj) = 1 / {ui — Ti + irj) , with 77 > 
infinitesimally small, and we hereafter set h = 1. This 
shorthand notation means, for example, that the (1,2) 
element of this matrix is: 

Gi,2(fc,c^) = (FM|c,,^_^G(a;)4^^_jFM) 

etc., so that all four possible combinations are considered 
at once. Using a Lehman representation,^^ it is clear that 



G\^y{k,uj) = 



(FM|c, 



\u){u\ci^,jFM) 



U! - 



Eu -I- i-q 



where T-L\u) = Eu\u) are the single-electron eigenstates 
in the sector of total z-axis spin NS — \. In other 
words, the poles of these quantities give all the eigenener- 
gies and the residues measure overlaps between the true 
eigenf unctions and the appropriate free-electron state 
i.J™),A = l,2. 

Our goal is to calculate exactly this Green's function 
matrix. We do this by using repeatedly Dyson's identity 
G{uj) = Go(w) + G{uj)VGq{uj), where G^{uj) = l/{uj - 
"Ho + ii]) ^-iid Ho is defined in Eq. ([9]) . 

Rotating to the diagonal basis cj: ^ ^ [see Eq. and 
back, it is straightforward to show that: 



Go(w)*t JFM) 



*t jFM)Go(^,w) 



where we introduce the 2x2 matrix: 



G^+^(fc,u;) G^ >ik,LO 



G'->{k,uj) Gi+\k,uj) J 
in terms of the known free-electron propagators: 



Gi^\k,u;) 



1 



± 



1 



(12) 



(13) 



(14) 



Using Dyson's identity once thus leads to the equation: 



G(k,uj)= 1 + (FMl^/j: ,G(uj)V^i\FM) Ga{k,uj 
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Since the objects appearing here are 2x2 matrices, the 
order of multiphcations in these equations is important. 
The action of V on ^1 |FM) is easily estimated. The 

xy exchange between spins has no contribution, since all 
spins are up. As a result, one finds contributions only 
due to the electron-spins exchange, resulting in: 



Gik,Lo) = Goik,Lo) 



Go{k,oj) (15) 



The first term in the square bracket comes from the di- 
agonal zz exchange and corresponds to a simple shift in 
the total electron energy (see below). The second term 
comes from the xy exchange, which allows one spin to be 
lowered by 1, while the electron-spin is flipped to a =1- 
This leads to a new 2x2 Green's function matrix: 

F{k,q,u) = Ee»«~^-^(FM|*,,^G(c^)*|__^5,^|FM>. 



Finally, the diagonal matrix: 



m 



cos ■ 







qya 



(16) 



appears in Eq. (|15p because the spin at site {ij) can be 
flipped by exchange with an electron present at either of 
the (jj, 1), (ij,2), (i - l,j, 1) or - 1,2) sites, which 
are displaced by ±ax/2 or ±ay/2 from it [see Eq. ^ for 
the definition of the phases in the Fourier transforms] . 
Given the structure of the Go matrix, it follows that 



JqS 



further simplify Eq. p5 



Ga{k,uj 
to: 



JoS) 



so we can 



G{k,uj) 



Jo \ ^ 

'n 2^ 



F(k,q,Lo)g{q) 



Go{k, Lu + JqS). 

(17) 



This shows that an equation for F is needed to solve the 
problem. Using Dyson's identity again, we find, after 
very similar kinds of manipulations, that its equation of 
motion is: 



is the one-magnon spectrum for this FM spin lattice. 

Eqs. (fT7|) and (fTS]) can now be solved to find G{k,uj) 
and F(k,uj). Note that usually there is an infinite se- 
quence of equations-of-motion connected to one another. 
In this problem, the series is truncated to just 2 equa- 
tions because of symmetries: since Sz^tot (which includes 
all spins and the electron) commutes with the Hamilto- 
nian, the evolution is always within the Hilbert sector 
with z-axis spin NS — 1/2. This only includes the states 
with all spins up and the electron with spin-down; or the 
electron has spin-up, and then one spin is lowered by one 
(one magnon is created in the system) . This also explains 
why generalizations to finite T (multiple magnons) or fi- 
nite electron concentrations are far from trivial: in those 
cases, more and more equations of motion are generated 
as the size of the relevant Hilbert subspace increases sub- 
stantially, and their solution becomes very difficult. 

Eqs. (fT7|) and (fT8|l can be solved analytically because 
of the simple structure of the g-matrix, which contains 
only trigonometric functions. In fact, if we also define: 



then we can "factorize" : 



sin2|ii 
sin^ 



m - q) = 9{Q)m + KQ)m 



(20) 



(21) 



smce cos 



(.Q-q)o. _ 



cos 



Qa 



COS ■ 



Qa • qa 

sm ^ sm 



We define the auxiliary quantities: 



(22) 



which is the only quantity we need to compute G, see 
Eq. (HT]), and 



/(fc,c.) = lE^(^'^'^)5(Q) 



(23) 



in terms of which we can rewrite Eq. (fT8|) as: 



F(fc, q, u) = 2 Jo^G(fc, Lo)g{q)Go{k - q,u; - r!g - JoS) 
- ^ E ^' - 9)Go(fc -q,u;- Q,- - JoS) (18) 



F{k,q,uj)^jQ -f{k,uj)g{q) ~ f{k,uj)g{q) 
+2SG{k,uj)g{q)] Go{k-q,uj-ng-- JoS) (24) 



where 



17,-=4J5fsin2^+sin2^" 
' V 2 2 



^ (19) 



Substituting this in Eqs. (|22|) and ([23)) . we obtain two 
linear (matrix) equations with unknowns / and /, and 
inhomogeneous terms proportional to G: 
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/(fc,u;)= -~Jof{k,uj) + 2JoSG{k,u) gii{k,uj) ~ Jof{k,uj)g2i{k,uj) 



and 



f{k,u))^ -Jofik.uj) + 2JoSG{k,uj) gi2{k,u)) ~ Jof{k,uj)g22ik,uj) 

r 



where we introduced the known 2x2 matrices: 
gii{k,uj) = -^^g{q)GQ{k-q,Lj-nif- JnS)g{q) (25) 

gi2{k,uj) = ^^g{q)Go{k-q,uj-n^~ JoS)g{cf) (26) 



In fact, lots of these matrices' elements are related to 
each other by various symmetries, so fewer than 16 actu- 
ally need to be calculated. Letting iV ^ oo, then each of 
these corresponds to a two-dimensional integral over the 
Brillouin zone. One integral can be performed analyti- 
cally, and the second we integrated numerically, therefore 
these matrices are easy to calculate. 



g2i{k,uj) = -^^g{q)Go{k-q,Lj-ni^- JoS)g{q) (27) 



92i{k,uj) 



N ^ 



These coupled equations are easy to solve and the re- 
sulting expression of / can now be used in Eq. (|17p to 
g{q)Go{k — q,LLj — f^g — JoS)g{q) (28) find the Green's function explicitly. The final result is: 



G{k,uj) 

where we introduced the matrix: 

M{k,uj) = l + Jo 



(Go{k,iu + JoS)'^ ^ -2S'Jo (l- (A?(fc,cj)' 



5ii(fc, w) - Jogi2(fc, w) ( 1 + Jag22{k,uj) ) g2i{k,uj) 

r 



(29) 



(30) 



Several comments are in order. First, the electron-spin 
exchange Jq appears in three places in Eq. (j29p . namely 
(i) in the overall shift by —JqS of the energy argument 
of the Gq^ term (first term in the denominator); (ii) as 
an overall factor for the "self-energy" (second term in the 
denominator) and (iii) as a shift by +JoS in the energy 
argument of the Gq functions appearing in the defini- 
tions of gii, etc. in the self-energy. The first and third of 
these are due to the zz exchange, which simply shifts the 
energy of the electron by ±JoS depending on whether 
its spin is parallel or antiparallel to the FM background. 
The second is due to the xy exchange, which facilitates 
the spin-flip of the electron. Therefore, the generaliza- 
tion to an anisotropic interaction is straightforward, for 
instance if Jo,_l — ^Jo,z then the self-energy is multiplied 
by A. The exchange J between spins appears only in the 
magnon dispersion Q^, which enters only the gn, etc. 
functions appearing in the self-energy. This is not sur- 



prising, since that self-energy term is due to contributions 
from one-magnon plus spin-up electron states. If the FM 
exchange between spins is longer range, one simply has 
to replace the magnon dispersion by the appropriate one. 

Eq. thus reveals that the free spin-down elec- 

tron state of bare energy e^j — JqS is mixed, through 
spin-flipping due to the xy term, with the continuum of 
one-magnon plus spin-up electron states of bare energies 



JqS 



fig-. It follows that if we are interested in 



having an infinitely long-lived quasiparticle state at low 
energies, then the electron-spin exchange must be anti- 
ferromagnetic Jq > (this is easily confirmed numeri- 
cally). We will focus on this situation in the following 
(the Jq < case is less interesting as the low-energy dy- 
namics is incoherent, with finite lifetime excitations). 

Let us now briefly comment on other generalizations. 
The reason we have equations for 2x2 matrices is that, 
in this problem, the electron unit cell has a two-site basis. 
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For a basis with n different sites per unit cell one would 
have similar equations but for n x n matrices, where all 
the electron hopping information would be encoded in 
the corresponding Go matrix. The dimensionality of the 
problem enters only in the sums over the Brillouin zone, 
i.e. in the munber of integrals to be performed. Slightly 
more complicated is the generalization of the electron- 
spin exchange to longer range. This leads to the ap- 
pearance of more auxiliary functions like / and /, i.e. 
one needs to solve a linear system with more unknowns 
in order to find the self-energy. The various auxiliary 
functions correspond to the different possible phase-shifts 
(analogs of the cos(qx^ya/2), sui(qx^ya/2) appearing in the 
g,g matrices). For example, assume that for the same 
system discussed here, the electron can visit all sites, in- 
cluding the spin sublattice. Also, assume that there is 
on-site exchange between the electron and the local spin, 
if the electron is on the spin sublattice, besides the ex- 
change discussed here. In that case, one has to deal with 
3x3 matrices (three-site basis) and there are three aux- 
iliary functions, two similar to the ones that appeared 
here, and one corresponding to zero phase-shift for the 
on-site interaction (because of the zero phase-shift, there 
is no auxiliary function proportional to sin(O)). 

We have checked explicitly that one can also include 
more complicated terms, for example electron hopping 
accompanied by a spin-flip coupled to a spin lowering of 
a nearby lattice spin, like ^ (such terms 

have important consequences, as we discuss below). In 
fact, we believe that essentially any problem from this 
class is solvable analytically, along these general lines. 



IV. RESULTS 

As already mentioned, we present results for the more 
interesting case of an AFM electron-spin exchange, Jq > 
0, whose low-energy state is an infinitely- lived spin- 
polaron. Also, the results are for ti = t2 = t (s-orbitals), 
and t will be used as the energy unit. As discussed previ- 
ously, for the simple nearest-neighbor hopping used here, 
the only difference for p-orbitals would be to shift the 
values of all momenta by -^(1,1). For holes, t changes 
sign with consequences discussed where appropriate. 

We begin by analyzing the dependence of the ground- 
state energy Eqs and quasiparticle weight on various pa- 
rameters of the problems. In panel (a) of Fig. [51 the full 
lines illustrate the dependence of Eqs on the electron- 
spin exchange Jo/t, as well as the spin value S, for a 
fixed spin-spin exchange J/t = 0.05. As expected, the 
GS energy is lowered as both Jq and S increase, mainly 
due to the favorable zz-exchange between the lattice of 
FM spins and the spin-down electron. 

The dashed lines in Fig. [D^a) show the asymptotic ex- 
pressions obtained by first order perturbation in the hop- 
ping t, in the strong-coupling limit Jo/t ^ oo. The agree- 
ment is very reasonable even for rather small Jq jt values, 
therefore it is useful to analyze this solution in some de- 
tail, to understand the nature of the spin-polaron. 

In the absence of hopping, i = 0, we can form two 
sets of translationally-invariant states of well-defined mo- 
mentum fc, which are ground-states of the electron-spin 
exchange Ticxc, namely: 



J 



1^,2)=^ 



1 



VN \ 4.S+1 



4Sc. 



|FM) 
|FM) 



(31) 
(32) 



For all of these states (A = 1,2), 

Hcxc|fc,A) =-Jo (•5+0 

The hopping T lifts this degeneracy. In fact, it is 
straightforward to show that (fc, A|T|fc', A) = while 

(fc, l|f|fc',2) = ^kk'^^H- Thus, within this level of 
perturbation theory, the eigenstates are 




and their corresponding energy is found to be: 
(v) 45 -I- i 3JS ( \\ 

The middle term is due to the spin-spin FM exchange 
energy lost because part of the wavefunction has one spin 
in the background lowered by 1. 

It follows that the perturbational prediction for the 
ground-state energy, shown in Fig. [21 is: 

In the limit of large J^jt the agreement with the ex- 
act solution is very good. In the limit Jq 0; and ig- 
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noring the small correction proportional to J (which is 
reasonable, given the small J/t value used), we see that 

> — 4i, whereas the true GS energy 



TP' 
tLi, 



(P) 



4S^ 



GS ^ ■^"45+1 

cannot go above —At (this is the free-electron GS en- 
ergy) . For large S the difference between these two values 
becomes negligible, however for S = ^ the difference is 
sizable, explaining the significant difference between the 
two curves in the Iow-Jq part of Fig. [2la). 

In Fig. Eljb) we show the exact GS energies for different 
values of J/t and S = ^. In the limit Jo/t ^ there is no 
interaction between the electron and the FM background, 
therefore the value of J is irrelevant. As expected from 
the discussion above, at large Jo/t values the GS energy 
increases linearly with increasing J. However, this is a 
small contribution for reasonably small values of J. In 
the following we fix J = 0.05t. 

In Fig. O^c) we plot the quasiparticle (qp) weights in 
the ground-state, defined as: 



Z 



AA = |(FM|c,,^o,A,ilGS)l' 



(36) 



Clearly, these quantities give the probability to find the 
electron with spin-down on the sublattice A = 1, 2, in the 
GS. As expected by symmetry, Zn = Z22 therefore only 
one is shown in Fig. [HJc), as a function of Jo/t. One 
immediate observation is that the qp weights saturate 
to finite values in the strong coupling limit Jo/t — s- 00, 
instead of becoming exponentially small, as is the case for 
typical polarons (where the electron binds phonons in its 
vicinity, due to electron-phonon coupling). The reason, 
of course, is that here the electron can bind a maximum 
of one magnon as opposed to an arbitrarily large number 
of phonons, as is the case with polarons for increasing 
electron-phonon coupling. 
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FIG. 2: (color online) (a) Spin-polaron ground-state energy, in 
units of t, as a function of Jo/t, for S = ^,1 and |,and J/t — 
0.05. The full lines show the exact results, while the dashed 
lines give the strong-coupling perturbational limit discussed 
in the text; (b) Spin-polaron ground-state energy, in units 
of t, as a function of Jo/t, for S = | and J/t = 0.05,0.25 
and 1.00; (c) Spectral weight Zn = |(FM|cgj JGS>p as a 
function of Jo/t for three values of 5* and J/t — 0.05. 



By direct comparison of Eqs. ([STll-dMl) and (g]), 
it is clear that in the strong coupling limit, the spin- 
down part of the low energy spin-polaron state is 
equal to \/4S'/(4S' + l)4.+,i|FM). In other words, for 

any momentum k of the spin-polaron, the probability to 
find the electron with spin-down and the FM background 
undisturbed is 45/ (45+1). By symmetry, it follows that 
Zu = Z22 — 25/(45 + 1) for all spin-polaron momenta, 
including the GS. This is in reasonable agreement with 
the values shown in Fig. [Ifc), given that even Jo/t = 7 
is not that large, therefore corrections beyond first order 
perturbation are not expected to be negligible. 

The spin-flipped contribution to the |fc, low-energy 
spin-polaron state also reveals very interesting physics. 
Focusing on the case 5 — 1/2, we can rewrite the spin- 
polaron eigenstates in the strong coupling limit as: 



(fc) 



|FM). 
.(37) 

The operator in the bracket is recognized as creating a 
singlet between the spin located at site ij and the electron 
occupying a state centered at site ij. This singlet prop- 
agates with momentum k through the FM background. 
The electron state that forms the singlet with the spin 

is found to be xh 



I t , t 

c _i . 1 ^ + e 2 c 



e~ 



* 2 c' 



t 



i.e. a superposition of 



-1 J,1,(T "T ^^"-1,2,(7! 

the four electronic sites surrounding the spin located at 
ij (site labeling is shown in Fig. [T]). For larger 5 the solu- 
tion is analogous, except one cannot speak of a "singlet" 
between a spin-i and a spin 5 > i. However, it can be 
shown that the entangled electron-spin state corresponds 
to a total spin 5—1/2 (not surprising, since this is total 
spin that minimizes the AFM exchange energy). 

In the ground-state, the coefficients are determined by 
the orbitals participating in the hopping. For our model 
we find for fc = that xj 



^i,j-l,2,(T- 



while for p-orbitals a similar calculation (or 



see k 

J 



J 



case) leads to xlj,^ ^ ~^^j,i,a " ^tj,2,a 
^'i-i j 1 (T + j-1 2 17- other words, the signs mirror the 
sign of the lob pointing towards the central spin site. 

If the charge carrier was a hole instead of an elec- 
tron, because t ^ —t the low-energy spin-polaron eigen- 
state is |A;, — ) (since now > 0). The eigenenergy 
is the same, but the orbital involved in the singlet is 



As a result 



for s-orbitals the GS orbital 



becomes xlj.„ = 4jA,a-4j,2^a+4-lJ,l_^a-4j~l,2,a^ ^'G- 

it has d-like symmetry. For p-orbitals it has p-like sym- 
metry again, but is orthogonal to the one listed above for 
the electron. 

In conclusion, the particular linear combination se- 
lected for forming the GS singlet (more generally, 5 — ^ 
state) with the central spin is determined both by the 
particular orbitals involved, and by the nature - hole or 
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electron - of the charge carrier. 

This solution is clearly analogous to the Zhang-Rice 
(ZR) singlel)^ but with some differences. For one, the 
singlet defining the spin-polaron is here propagating in a 
FM, not AFM background. Second, while the ZR singlet 
involves a d-wave like linear combination of electronic 
orbitals, here the combination depends on the details of 
the model considered, as just discussed. 

The bigger difference, of course, is that here we have 
spins at the Mn-like sites. Of course, spins arise from hav- 
ing some atomic orbital partially filled, and one can talk 
about a well defined spin when the number of electrons 
(holes) in this shell cannot change. As noted above, we 
can easily generalize our model to allow the extra electron 
(hole) to hop onto the Mn-like sites, adding a Hubbard- 
U penalty and/or Hundt's exchange as well, if desired. 
What we cannot do, at least so far, is solve exactly the 
more general model where electrons (holes) that are cur- 
rently locked onto the Mn-like sites and constitute their 
spins, are allowed to hop to the other sublattice, so that 
a spin less than S is left behind. The difficulty is simple 
to see: even if one adds a charge-gap A, i.e. an energy 
penalty to move electrons (holes) from the Mn-like sites 
to the 0-like sites, in any eigenstate there will be some 
finite probability to find any number of O-like sites occu- 
pied and the wavefunctions become too complicated. To 
be more precise, one can still find easily the equivalent 
of the FM-background state in this case, i.e. the ground- 
state in the Hilbert sector of z-axis spin NS. The case of 
spin NS + n^, when any number n of electrons (holes) 
with spin parallel to the background have been added in, 
is also trivially solved (because Hubbard on-site repul- 
sion does not act between parallel spins, therefore both 
cases are essentially without interactions). However the 
problem corresponding to NS — 1/2, i.e. for adding a 
spin-down electron (hole) to the FM "background", in 
other words the equivalent of the simpler problem inves- 
tigated here, has proved too complicated for us so far. 

However, even the asymptotic limit of our simplified 
model still provides a very important insight, namely 
that the phases of the electronic orbitals locked in the 
singlet with the central spin vary inside the Brillouin 
zone. In the ZR model& these phases are assumed to be 
locked to their GS, d-wave symmetry values irrespective 
of the momentum of the ZR singlet. As is well known, 
that leads to problems with normalization of the result- 
ing states in some regions of the BZ. What our simpler 
but exact solution reveals is that this is not correct: the 
phases of the orbitals involved in the singlet vary at dif- 
ferent fc-points. This insight might help improve the de- 
scription of the ZR singlet away from the fc = region. 

In terms of fc-dependent properties, of course we can 
extract not only the GS state, but the entire dispersion 
of the spin-polaron, by focusing on the lowest (discrete) 
eigenstate of momentum k. In standard polaron physics, 
the polaron bandwidth is expected to become smaller 
(corresponding to larger effective polaron mass) as the 
electron-phonon coupling increases and more phonons are 




(k ,k ) 

X y 

FIG. 3: (color online) Spin-polaron dispersion along lines of 
high symmetry in the Brillouin zone. The results (symbols) 
correspond to 5 = ^ and J/t = 0.05 and, from top to bottom, 
Jo/t = 1,2,3,5 and 7. The top two cases have been plot- 
ted separately to avoid overlaps. The lines show the strong- 
coupling prediction -EI''^ of Eq. I|34[) . 

tied into the polaron cloud. Here we do not expect to 
have this problem since a maximum of one magnon can 
be bound to the electron, as already discussed. In fact, 
Eq. ([M)) reveals that the spin-polaron mass becomes, in 
the strong-coupling asymptotic limit: 



where m is the bare electron mass. This is interesting as 
it is determined only by the value S of the spins in the 
FM background, and independent of the coupling. The 
largest increase of 6/5 = 1.2 is for spins S = 1/2, showing 
that the spin-polaron remains very light. This is not 
surprising, given its propagation in an FM background. 

In Fig. [3] we plot the spin-polaron dispersion (lowest 
eigenstate of momentum k) along lines of high symmetry 
inside the BZ. The full lines show the strong-coupling pre- 
diction E'^y' of Eq. ([M]) - again, we assume s-orbitals 

and an electron as the charge carrier. The agreement 
is very good for larger J^/t values, suggesting that the 
paradigm of the spin polaron as a mobile singlet of the 
electron in a special state and the spin at its center, must 
be a good description even for moderate values of Jo/i. 
Even for J^/t = 1 the agreement is quite reasonable near 
the center of the BZ, but not at high k values. This dis- 
crepancy can be understood easily as well. As discussed, 
the spin-polaron states come from hybridizing the non- 
interacting spin-down electron states of energy e^: — JqS 
with the spin-up electron -I- one magnon continuum of 
energies ej:_^ + JqS + Jig. The hybridization will push 
the discrete state to lower energies, giving rise to the 
spin-polaron band, however it cannot change the loca- 
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tion of the continuum. If we ignore the fig- term, given the 
small J/i, then we find that this continuum must start 
at — 4t + JqS (inclusion of adds a small, fc-dependent 
correction to this value of the continuum band-edge, see 
results below). The spin-polaron discrete state cannot 
overlap with the continuum, therefore when Jq is small, 
one expects the polaron dispersion to flatten out just be- 
low the continuum band-edge, which is precisely what we 
see for J^/t — \. As Jo/t increases the continuum moves 
to higher energies while the spin-polaron band moves to 
lower energies and becomes fully visible. 

The flattening-out of the polaron dispersion just un- 
der the continuum is typical polaron physics at weak 
coupling. Of course, for conventional polarons the con- 
tinuum starts at Q above the polaron ground-state, where 
r2 is the energy of the Einstein phonons. This limits the 
polaron bandwidth, at weak coupling, to be precisely $1. 
For the spin-polaron studied here, the bandwidth has lit- 
tle to do with the magnon energy (the FM magnons are, 
in fact, gapless). Instead, the relevant energy scale comes 
from the zz component of the exchange energy between 
the electron and the FM background, JqS. 

Another typical expectation for polaron physics at 
weak couplings is that of very small spectral weight in the 
BZ regions where the polaron dispersion lies just below 
the continuum. This is because here the largest contribu- 
tion to the polaron comes from continuum states where 
the electron is spin-up. Consequently, the probability to 
find a spin-down electron, measured by the qp weights, 
becomes very small. This is indeed confirmed in Figs. 21 
where we show the spectral weight 



Axx{k,uj) 



1 



(39) 



along the = ky cut in the BZ. Again, by symme- 
try we expect An — A22- The top panel shows the 
spectral weight on a linear scale. The low-fc part of 
the spin-polaron band is clearly visible, however the 
region where it flattens just under the continuum has 
very little weight, and is not visible on this scale. In 
order to make low-weight features near the continuum 
more visible, instead of Aii(fc,a;), in the bottom panel 
we plot tanh[Aii(fc,a;)/0.3]. As a result, all regions 
with weight An > 0.7 are mapped into black, whereas 
weights below this follow a fairly linear scale down to 
white. This explains why the spin-polaron peak seems so 
wide now, although in reality it is a Lorentzian of width 
r] 0. The dashed line shows the continuum band-edge, 
minq-[e^_^ -1- JqS + fi^]. As expected, the spin-polaron 
dispersion flattens out just underneath it. Because its 
weight decreases so fast as k increases, it is impossible 
to see it even on this scale for larger k. The contin- 
uum above is also more clearly seen for smaller values 
of /c, with most weight around eg — JqS where the free 
spin-down electron would have its bare energy. The lit- 
tle weight that seems to "seep" below the dashed line at 
low-fc is due to the finite ry used, we have checked that 




(0,0) 



(k,k) 




(0,0) 



(kx ,ky) 



(7t,Jl) 



FIG. 4: Top panel: Contour plot of Aii{kx,ky,u!) along the 
kx = ky cut in the Brillouin zone. Only the spin-polaron 
state is visible on this scale. Bottom panel: same as above, 
but plotted as tanh[^ii (fca;, fcj/, a;)/0.3] so as to make low- 
weight features visible. The dashed line shows the expected 
location of the lower edge of the continuum. Parameters are 
Jo=t = l,J = 0.05, S = 0.5 and r? = 0.02. 



indeed the continuum appears at the expected value. 

As Jo increases and the spin-polaron band moves well 
below the continuum, we expect to see large qp weight for 
the spin-polaron at all k since the strong-coupling limit 
predicts a fc- independent qp weight. This is indeed the 
case, as shown in Fig. [5] for Jo/t = 5. Note that the 
spectral weight is now shown on a linear scale. The low- 
energy polaron band has the same dispersion as shown in 
Fig. [21 and indeed it has a fairly constant qp weight. One 
initial surprise is to see another coherent, infinite lifetime 
state just above it (and in fact degenerate with it along 
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FIG. 5: Contour plot of A-i_i{kx,ky,uj) along several cuts in 
the Brillouin zone. The dashed line shows the expected lo- 
cation of the lower edge of the continuum. Parameters are 
t = 1, Jo = 5, J = 0.05, 5" = 0.5 and r/ = 0.1. 



the (tTjTt) — (tTjO) cut, leading to a doubled qp weight 
there), however this is to be expected. The spin-polaron 
also lives on a lattice with a two-site basis, so one expects 
to see two bands for it, if both happen to fit below the 
continuum, as is the case here and for larger Jq values. 
This is also in agreement with the strong-coupling limit 
which also predicts the two bands, see Eq. From 
there, we see that this higher energy spin-polaron state 
can also be thought of as a propagating singlet, where the 
singlet again involves the electron in a particular wave- 
function and the spin at its center. The wavefunction 
involves the four sites surrounding the spin, but with dif- 
ferent fc-dependent phases (or symmetry) than for the 
low-energy branch. 

Besides these two spin-polaron bands. Fig. \5\ reveals 
another dispersing feature at much higher energies, of 
order —2t. That this is also a discrete (infinitely lived) 
state is demonstrated in Fig. [HJa), where we verify its 
scaling with rj: doubling rj halves the height of the peak 
and doubles its width. Thus, this is a Lorentzian of width 
rj, and becomes a delta function as 77 ^ 0. Note that the 
continuum is independent of rj and begins at the expected 
value, indicated by the vertical line (the inset focuses 
on this feature and demonstrates this agreement more 
clearly). The continuum has very low weight, explaining 
why it is mostly invisible on the scale of Fig. [51 

Interestingly, we only see this higher-energy discrete 
state in some regions of the BZ - for example it is absent 
along the (0, 0) — (0, tt) hne (to be more precise, we have 
searched with an rj as low as 10~^ and did not see any 
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FIG. 6: (color onhne) (a) Spectral weight Aii{^, ^,uj) as 
a function of energy, for rj = 0.01 (dashed line) and rj = 
0.02 (fuU line), (b) ^ii(7r, f,a;) (full line) and Aii{^,n,cj) 
(dashed line) vs. cu, for rj — 0.02. In both cases the vertical 
lines show the expected location of the lower band-edge of 
the continuum. The insets show the same quantities, but 
with the focus on the lower band-edge of the continuum. The 
parameters are t = 1, Jo = 5, J — 0.05 and S = i. 



features just below the gap. Of course, this does not rule 
out a state with an extremely low qp weight). The ap- 
parent "disappearance" near (tt, tt) is due to vanishing qp 
weight, but the state exists in that region. Another inter- 
esting feature is demonstrated by Fig. [Hfb), which shows 
Aii{k, Lu) at two points that one would normally consider 
equivalent, namely (tt, ^) and (f ,7r). While the former 
shows a big qp peak just above —2t, the later shows no 
weight at this energy. The explanation is that by sym- 
metry one needs to have Aii{kx,ky,u!) = A22{ky,kx,uj), 
since "1" and "2" refer to states on the x, respectively 
y-rungs. However, for kx ky, there is no requirement 
that Aii(kx,ky,u!) = Aii{ky,kx,uj) and indeed we see 
that this is not the case. Fig. [Bfb) thus suggests that the 
disappearance of this qp state along some directions is 
due to symmetries, which result in orthogonality between 
the spin-down free electron state and the true eigenstate. 
Polarization-dependent spectroscopies should be able to 
detect these variations between An and ^22- 

Regarding the origin of this high-energy qp state, this 
seems to be standard polaron physics. Holstein polarons 
at strong couplings are also known to have a so-called 
"second bound state".— It is essentially an internal ex- 
cited state of the polaron, possible if the interaction and 
so the binding energy with its cloud is strong enough. 

These two examples, together, are enough to give us 
intuition about the evolution of the spectral weight in 
such models as Jq increases. For small values only the 
lowest spin-polaron band is visible in regions where it 
is well below the continuum. In the regions where the 
continuum forces it to flatten, the qp is extremely small. 
As Jo increases, more and more of the low-energy spin- 
polaron band emerges below the continuum. Then the 
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FIG. 7: (a) Illustration of exchange between the charge carrier 
and the impurity spin (depicted as a box), arising from virtual 
hoping of the charge carrier to the spin site, and back, (b) 
The same type of process can result in charge carrier hopping 
with spin flip, if the charge carrier moves to another site. 



upper spin-polaron band starts to emerge as well and 
the part that is below the continuum is clearly visible, 
the rest is flattened and with low qp weight. This is 
the case e.g. for Jq = 3 (not shown) where the entire 
lower band and about one third of the upper band are 
visible below the continuum. For even larger Jq both 
spin-polaron bands are visible, and at yet higher values, 
more bound states can split-off from the continuum. The 
number of qp bands and their weights is therefore very 
sensitive to strength of the coupling. 

Before concluding, we must make one more important 
comment. In the model considered here we assumed that 
the charge carrier interacts with the spins only through 
exchange. This exchange comes from virtual hopping of 
the charge carrier to the spin site and back. Fig. [7[a) 
illustrates one such possible scenario, leading to AFM 
exchange if a single half-filled orbital gives rise to the 
impurity spin. However, the charge carrier could hop off 
to any of the nearest neighbor sites of the impurity spin, 
since it does not necessarily have to return to the original 
site. Such processes are illustrated in Fig. (Tljb) and result 
in effective charge carrier hopping with possible spin-flip 
correlated with a spin flip of the impurity spin. The 
energy scale for such processes is the same as for the 
AFM exchange, however their sign can be either negative 
or positive, depending on whether the hopping integrals 
along the relevant links have the same, or opposite signs. 
This, of course, depends on the orbitals involved both at 
the charge carrier and the spin site. 

Such hopping + spin-flip processes can be treated ex- 
actly in our approach, as already stated. They can have 
dramatic effects, as we illustrate now with a simple ex- 
ample. Consider the spin S at site (ij) and its four 
neighboring charge sites, which we label anticlockwise as 
1 = 2 = (zj,2), 3 = (i-l,j,l) and4 = {i,j-l,2). 

Four symmetric linear combinations can be built from 

these, namely Xsa- = \{cia + C2a + C^a + Cia)] Xda = 
i(cicr-C2<T + C3CT-C4cr); Xpla = ^ (cia + C2a - - Cia) and 

Xp2a = ^(cicr " C2a - c^a + C4a) ■ Cousidcr the impurity 
spin lowering part of the Hamiltonian we have worked 
with so far, JqS^ {c\^cii + c\^C2i -\- c\-^C3[ + c\^cai) — 



.e. all these 

symmetrized orbitals interact equally strongly with the 
impurity, and could participate in the "singlet". Which 
happens to be the GS is chosen by the signs of the hop- 
pings and charge carrier, as discussed. If hopping + spin 
flip processes are included, the results are very different. 
Assume for example that all hopping integrals between 
the 4 sites and the impurity site have the same sign. In 
this case, JqS^ c\^cii JqS^ {c\^ + 0^2] + + C4-f)cix, 
since as illustrated in Fig. [Tj^b), the charge carrier re- 
moved from site 1 can emerge at any of the spin's neigh- 
bor sites. The same is true for all other terms, and adding 
them together, we obtain a total equal to ^JoS~^^^Xsi- 
In other words, in this case only the s-symmetry orbital 
can participate in the singlet (which state is selected de- 
pends, obviously, on the relative signs for the hopping 
between charge sites and the impurity site). As a result, 
we expect to see a single spin-polaron state at fc = 0, cor- 
responding to a "singlet" involving this allowed orbital, 
and the other spin-polaron states discussed for the simple 
model will vanish from the spectrum. 



V. SUMMARY AND CONCLUSIONS 

To summarize, we claim that essentially any problem 
that involves a charge carrier interacting with a fully po- 
larized FM background can be solved exactly, irrespective 
of the complexity of the sublattices involved or the range 
of the exchanges or other complications. Of course, cal- 
culations become more complex as one makes the model 
more complicated, but the exact solution following the 
approach suggested here should always be possible. 

We exemplified it on a relatively simple, yet interesting 
two-dimensional case which we believe illustrates most of 
the physics that can be expected in such systems. It al- 
lows us to elucidate the nature of the spin-polaron quasi- 
particle, which at strong coupling is a singlet (more gen- 
erally, the maximally polarized state with spin — i) 
similar to the Zhang-Rice solution. However, we show 
that the orbitals participating to the electronic wavefunc- 
tion that locks into the singlet with the lattice spin, have 
phases that vary as a function of the singlet's momen- 
tum. In the ground state they have the signs consistent 
with the symmetry of the orbitals involved and also de- 
termined by the type of charge carrier. Moreover, since 
on complex lattices one can form multiple such linear 
combinations which can participate in the singlet, it is 
possible to see more than one spin-polaron band below 
the continuum. If however hopping -I- spin-flip terms are 
included, some of these spin-polaron bands will be re- 
moved because of symmetry considerations. 

We hope that these insights may help solve some of 
the known issues regarding the normalization of the ZR 
singlet, which is an essential ingredient in the cuprate 
physics. Extensions of the solution to AFM backgrounds 
are very desirable, although clearly much more compli- 
cated. We plan to investigate such problems next. 
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As far as uses of doped FM insulators for spintron- 
ics applications are concerned, this work suggests that 
materials with an AFM exchange Jq between the charge 
carriers and the local spins could be very interesting. For 
instance, consider a pulse of unpolarized charge carriers 
created in such a medium, for instance by optical means. 
If they are placed in an electric field, the charge carriers 
with spin parallel to the FM background will propagate 
with a different speed than the charge carriers with spin 
antiparallel to the FM background, since the latter are 
dressed by magnons and become heavier spin-polarons 
while the former propagate as bare particles. Both types 
propagate coherently (at least at T = 0) so this suggests 
that in time a pulse of unpolarized carriers will separate 
spatially into two spin-polarized pulses traveling at differ- 
ent speeds. Ingenious use of such difference between the 
two spin polarizations may open the way towards using 



such materials as sources or detectors of spin-polarized 
currents, which are essential components for spintronic 
devices. Of course, issues such as the effect of finite-T 
and finite charge carrier concentrations need to be un- 
derstood first, however this seems to be a promising line 
of investigation. 
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